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We consider inflationary models in which vector fields are responsible for part or eventu- 
ally all of the primordial curvature perturbation £. Such models are phenomenologically 
interesting since they naturally introduce anisotropics in the probability distribution 
function of the primordial fluctuations that can leave a measurable imprint in the cosmic 
microwave background. Assuming that non-Gaussianity is generated due to the super- 
horizon evolution, we use the SN formalism to do a complete tree level calculation of 
the non-Gaussianity parameters /nl and tnl in the presence of vector fields. We isolate 
the isotropic pieces of the non-Gaussianity parameters, which anyway have contributions 
from the vector fields, and show that they obey the Suyama-Yamaguchi consistency re- 
lation rj|£ ^ (f /nl) 2 ' Other ways of defining the non-Gaussianity parameters, which 
could be observationally relevant, are stated and the respective Suyama-Yamaguchi-likc 
consistency relations are obtained. 

Keywords: Non-gaussianity; statistical anisotropy; vector field models. 
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1. Introduction 

The study of non-Gaussianities in the primordial curvature perturbation £ is a 
subject of major interest in modern cosmology because the evaluation of the non- 
Gaussianity (NG) parameters provide criteria to discriminate among the many mod- 
els proposed to explain the origin of the large-scale structure that we observe today^. 
Particular attention has be en p aid to the study of a consistency relation between 
the /nl and tnl parameters^!, the Suyama-Yamaguchi (SY) consistency relation 
(tnl ^ (|/nl) 2 at least at tree level), because its violation would rule out many 
of the popular models of inflation. So far, it has been shown that this consistency 
relation works, in principle, on models that only include scalar fields under just a 
few assumptions, but there is no any conclusive result for models that include other 
type of fields, for instance, vector fields Vector field models are particularly 

interesting because they are suitable candidates to explain the apparent violation of 
statistical isotropy observed in recent analysis of data from the WMAP satellite^, 
given that VF define inherently a preferred dire ction f or the expansion, for the dis- 
tribution of primordial fluctuations, or for botr pEE Because of this reason, it is 
pertinent to study consistency relations in models that include scalar fields as well 
as vector fields as these not only include NG parameters but also the amount of 
primordial statistical anisotropy <^K 

The main purpose of this letter is to find out how the well known SY consistency 
relation, valid for models involving only scalar fields under just a few assumptions, 
is modified when including VF. To this end, we do a complete tree level calculation 
of the /nl and tnl parameters based on the SN formalism^ and concentrate on 
the isotropic pieces of these parameters since these are the ones that are currently 
constrained by observations. Under pretty the same few assumptions made in the 
scalar fields case, and although the VF do contribute to the isotropic pieces of 
the NG parameters, we prove that the SY consistency relation is still obeyed, at 
least at tree level. Finally, the modified SY consistency relations for other ways of 
defining the NG parameters, which could be observationally relevant when taking 
into account the level of statistical anisotropy, are obtained. 



2. The curvature perturbation 

In the presence of VF, part or even all of the primordial curvature perturbation 
£ can be generated by vector field perturbations. In such scenarios, the curvature 
perturbation can be calculated through the SN formalism as an expansion in the 
perturbations of the fields a few e-folds after horizon exitK 

C = SN = N a SA a + -N ab SA a 5A b + -N abc 8A a 6A b 5A c + • • • , (1) 

a VF naturally generate scale-dependence and, therefore, the proofs in Ref. El do not apply to those 
scenarios where VF contribute to the generation of the primordial curvature perturbation. 
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where N is the number of e-folds during inflation and N a b... c are the derivatives of 
N with respect to A a . In the previous expression, and in order to avoid notational 
clustering, we have assumed that there is only one scalar and one vector field; thus, 
SA a = 6cf) for a = 0, and 8A a = SAi, the latter being the spatial components of the 
vector field, for a = i = 1, 2, 3^1 where these field perturbations are calculated in a 
flat slicing (the threading must be the comoving one). The extension of the results 
found in this letter to the case of several scalar and vector fields is straightforward. 



3. Calculation of the correlators 
3.1. The Power Spectrum 

The power spectrum (PS) of the primordial curvature perturbation £ is defined as: 
(C(fci)C(fe)) = (2^) 3 <5(fc 12 )P c (fc 1 ) (2) 

= (27r) 3 <5(fc 12 )|f7' c (fc 1 ), (3) 

where 7\(fc) is the dimensionless power spectrum, and fci 2 = fci + fc 2 . The statistical 
anisotropy in £, due to the existence of a preferred direction h, can be parametrized 
by a function as follows^ 

P c (k)=P™(k)(l+g c (k-h) 2 ), (4) 

where P^ s °(k) denotes the isotropic part of the PS and the vectors k and n are of 
unit norm. 

The first assumption we will make is that the NG in the field perturbations is 
negligible, i.e., the NG is generated due to the superhorizon evolution. The second 
assumption is to consider that the field perturbation^] are statistically isotropic, 
which is equivalent to the fact that there is no correlation between scalar and 
vector fields and that the expansion is isotropicP. This may be reasonable as a good 
approximation in several caseJ^, although it must be proven in each particular case, 
and is essential in the development of this letter because, otherwise, there would 
not exist any isotropic contributions to the NG parameters which are what really 
is constrained by observations. We shall adopt the notation 

(SM^SAbfa)) = (2n) 3 S(k 12 )U ab (k 1 ) , (5) 

with the components 

(Mo(fci)M (fe)) = (27r) 3 5(k 12 )P s 4k 1 ) , (6) 
(8Ai$i)SMh)) = (27r) 3 «5(fc 12 )n iJ (fc 1 ) , (7) 
and zero otherwise. In the latter expression: 

ny(fc) = u^(k)p + (k) + n°? d (k)p_(k) + n:° ng (fc)p long (fc) , (8) 



b We refer, in the case of VF, to the scalar perturbations that multiply the respective polarization 
vectors in a polarization mode expansion. 
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where II^ cn (fc) = 6,, /,-,/,■,. II° dd (fc) = ie ijk k k , and II*° ns (£) = hkj, P long (k) being 
the longitudinal power spectrum and P+ (k) and P_ (k) being the parity conserving 
and violating power spectra respectively^! The fact that none of Pg^k), P+(fc), 
P_(fc), and Pi ons (k) depend on the direction of the wavevector is a consequence of 
the second assumption. 

At tree level, the PS isP^ 

P c (k) =N a N b n ab (k). (9) 

For concreteness, we shall restrict to the case of vector field perturbations which 
preserve parity, then, we just keep the P+(k) and P\ ng(k) parts of the vector field 
power spectra. We also suppose that all the <L4j perturbations evolve in the same 
way after horizon crossing so that they have the same spectral index in their power 
spectra. In such conditions, the longitudinal and the parity conserving power spectra 
are related by P\ on g(k) = rP + {k) where r is a scale independent number. With these 
considerations in mind, the vector power spectra are: 

Htf (*) = [<% + (»• - p +( k ) ■ (10) 

Our next steps are to evaluate the PS by employing the SN expansion and to 
separate the isotropic and anisotropic pieces obtaining in this way: 

P c (k) = N a N b P ab (k) + (r- 1) (Nik) 2 ' P+(k) = Pr(k) (l + 9c (h t h) 2 ) , (11) 
where 

Pf°(fc) = N a N b P ab (k) , (12) 

P c aniso (fc) = (r- 1) (nA) 2 P+(k) , (13) 

hi = Nil ' {NjNj) 1 / 2 and the relation between the (r — 1) factor and the anisotropy 
parameter in the power spectrum is given by: 

gC = (r-l) W+ . (14) 
yC 1 'N a N b V ab 1 ' 

In the latter expression, the dimensionless power spectra V ab assume the values 
Poo = Vs4> and Vij = S l] V + . 

3.2. The Bispectrum. 

The bispectrum (BS) of the primordial curvature perturbation is defined as: 

3 

([[ah)} = (27r) 3 S(k 123 )B ( (k ll k 2 , k 3 ) (15) 
»=i 

47T 4 - w _ 

= (2^f5{k 125 )-^B c {k u fc 2 , k 3 ) . (16) 



November 26, 2012 1:14 WSPC/INSTRUCTION FILE ws-mpla 



The Suyama- Yamaguchi consistency relation in the presence of vector fields 5 

Using the SN expansion we get the tree level BS: 



B c = N a N b N cd n ac (A:i)n M (fc 2 ) + 2 perm. 



(17) 



Expanding the above by employing the expression in (| 10[) and separating the tree 
level BS into the isotropic and the anisotropic pieces we get 



Bf°{ki, k 2 , fc 3 ) = NaNbNdVacVbd ^(^m) -3 , (18) 
B% niso (k!, ha, k 3 ) = (r - l)N a N l N bj V ab T+ (kik m )- 3 (k { i )t k (l)j + k (m)l k (m)3 ) 
+ (r- l) 2 NiN k N jn Vl ^(kik m )-% )i k Wj k {m)k k (:m)n • (19) 



3.3. The Trispectrum 

Finally, we evaluate the trispectrum (TS) of the primordial curvature perturbation: 



(U_C(h)) = (2n) 3 S(k 1234 )T ( (k 1 , k 2 , h, k 4 ) 



(27T) 3 J(fci234) 7 3 1 3 7 3 Tt(kl, k 2 , fa, fc 4 ) 
">1 "> 2 ^23 



(20) 



(21) 



Once again, we use the SN expansion to get the tree level contribution to the TS, 
which results in 



= N a N b N cd N ef 



Uac(ki)U be (k2)Tldf(ki3) + 11 perm. 



N a N b N c N def n ad (/c 1 )n he (/c 2 )n c/ (A:3) + 3 perm. 



(22) 



(23) 



As we are interested in the tnl amplitude, we only need the term T™ L given that 
the term Jj? NL enters in the definition of the ^nl parameter. The decomposition of 
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T TNL in isotropic and anisotropic pieces is given by 
Tl s ° imL) = N a N b N cd N ef PacV be r df (hkmk ls )- 3 , (24) 

Z<m.s^m,Z 

T a„iso(r NL ) = (r _ 1)N . NbNjdNefVbeV(JfV+ (klk m k ls )- 3 Ck {l)t k (l)j + k {m)i k {m)j ) 

l<Cm,s^m,l 

+ (r - l)N b N d N el N fj V be V df V+ Y, ( k lkmk ls r 3 k {ls) ik {ls)j 

l<Cm,s^m,l 

+ {r-l) 2 N l N k N ]d N lf V df Vl Y. (kik m k ls )- 3 k {l)l k (l)j k 

+ (r - l) 2 NiN d Nj k Nf n VdfV% ^ (hk m ki s )~ k^kky^ik^ik^j + 

+ (r- lfN % N k N om N hn Vl ( k lk m k ls y 3 k {l)l k(i)jk^ m )j~k^ m ^f l k(i s ^ m kfi s y n . (25) 



4. The NG parameters /nl and tnl 

In this work we consider the definition of the non-gaussianity parameters /nl and 
tnl as follows 

B c = ~fMPdh) P dh) + 2 perm.) , (26) 
o 

T c = r NL ( J Pc(fci)^c(fc2)Pc(fci4) + 11 perm.) . (27) 



Notice that, in the definitions above, we have used the full PS including both scalar 
and vector field perturbations which, accordingly, include isotropic and anisotropic 
contributions. Nevertheless, it is usual to find in the literature a definition which 
considers only the isotropic piece of the P£P, so the NG parameters are given by 

B ( = !/nl(^c so (£i)^c S °^ 2 ) + 2 P erm -) and T C = ^nl(^ so (^i)^c SO (^) P C S °(^ 4 ) + 
11 perm.). We can relate both definitions using the expression in (fTTj) for the PS: 

•^ NL = i I Vl = i 7 v-2 r~~ /nl , ( 28 ) 

1 + Xi l + £ i= i«M5c) 

t nl = t— — T NL = — -=g J nl , (29) 

1 + X2 1 + £,•=!** (fife)* 



where the parameters rrii and li are obtained from the expansion of the products of 
the PS in (fTTj) which are in the definitions of the NG parameters. They depend on 
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the derivatives of N and the momenta ki. Explicitly: 

P f (fc 1 )P c (fc 2 ) + 2perm. = + [pf (ki)P£ (.k 2 ) + 2 perm.] 

„ „ 12l<m( k l km ) 3 (fc(i)jfe(;) ;/ - + fe( m )jfc( m )j) 

1 + 9C n ^ 12i <m (hk m )-z 



2 - * « „ Y,l <m ( k l km ) k (l)i k (l)j k (m)k k (m)r. 

+ g Q n %nj n k n n — /, , \ , — 

^ E;< m (kk m ) 3 



(iV a iV b P ab ) 2 ^(fc z fc m )- 3 , 

(30) 



P c (fc 1 )P c (fc 2 )P c (fc 13 ) + 11 perm. = (1 + X2 ) [pf^fc^Pf "(fcji^fcia) + Uperm 

„ „ E/< m s^m i ( k l k rn.ki s ) 3 { k (l)i k (l)j + k (m)i k (m)j + k (ls)i k (ls)j) 

1 + g c n 4 nj ■ = rfcjfc fc , w 3 

l^l< mi s^m.l\ K l KmK ls> 

2 E/< m s^m li k l k m k ls) 3 &(;)jfe(i)j&( m )fcfe( m )n + k (ls)i k {ls)j ( k (m)k k (m)n + ^(i)fe^(/)n) 
+ ,g c n^hkhn ■ — - 

3 Y^l<m,s=£mA k l k m k ls) 3 ' k (l)i k (l)j k (m)k k (m)h k (ls)t k (ls)n 

+ g ( n inj n k n h ntn n ■ = 3 

ANaN b T ab f J2 ( k lkrnk ls )- 3 . (31) 

Certainly, both definitions of the NG parameters are approximately equal if the 
corrections arising from the existence of some statistical anisotropics were negligible. 
Otherwise, such corrections would be very important from the observational point 
of view. 

4.1. The /nl parameter 

For the /nl parameter we can use the results (|T2|) . (fl~3|) , (IT8|) . and (p~9|) to sepa- 
rate the isotropic and the anisotropic contributions as follows: /nl = /nl + /nl = 
B'' c so /(P c P c + 2p.) + J B^ niso /(P c P c + 2p.). The labels I and II are related to the 
isotropic and the anisotropic pieces of the bispectrum in (fT8|) and (|T9|) respectively. 
Using this decomposition, we can compare the / NL and /n l contributions by eval- 
uating its ratio 



/nl 



1=1 



NaNbjN.N.NdVabPcd El<m (Mm)' 3 (*(J)i*(J)j + k 

(m)ifc(rn)j ) 

= 5C N a N b N cd N k N k P ac P bd Ei<m( fc i fc »0" 3 
a (N a N b P ab ) 2 NiNkN jn T,i< m ( k i k m)' 3 k { i^ k (i)j k 

9C (NiNi) 2 N a N b N c dPacPbd E^JMnTF 3 ' 

In the latter expression, we have used the expression in (fT4")l to write the series 
expansion in terms of the level of statistical anisotropy g^. The coefficients / can 
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be read off directly from the latter expression. Using the ratio /nl//nl> we can 
write the total / NL as: / NL - (1 + £i)/ NL = (1 + + Xife where is 

the piece of /nl that depends only on the isotropic BS in (IT51) and the isotropic 
PS in (|12[) . A very appealing feature that we can see from the previous analysis 
is that the ratio of the terms related to the anisotropic and the isotropic parts of 
the NG parameter /nl depends on the parameters of the model; for instance, it 
depends on the value of the ratio r, the level of statistical anisotropy g^, and the 
shape of the triangle defined by the ki momenta. Moreover, interestingly enough, 
for some values of and for some configurations of the momenta, the anisotropic 
contribution could be of the same order of magnitude as the isotropic piece and 
cannot be neglected. 

4.2. The tnl parameter. 

For the tnl parameter, we apply the same decomposition related to the isotropic 
and anisotropic pieces that we did before with the /nl parameter: tnl = t nl+ t nl = 
Tl so /(P ( P ( P c + lip.) + T* niso /(P c P c P c + lip.). We follow the same steps that we 
did in the case of the /nl parameter and write the ratio of both contributions as a 
series in powers of the g^ parameter: 

NL (=1 

m)i^(m)j ) 

' ^ 



+9C 



, 2 

+9( 



{NiNi)NaN b N cd N ef Vacr be V d f J2i <m , S7 t m ,i(kik m ki s )- 3 
(N a N b V ab )N b N d N el N fj P be V df T,i <m , s ^ m ,i(hk m ki s )- 3 k {ls)i k ils)j 
(N i N i )N a N b N cd N ef VacVbe'Pdf T,i <m , S7 t m ,i(hk m k ls )- 3 

(N a N b T ab ) 2 NiN k N jd N nf V df EKm, S #m,;( fc i fc " fc is)~ 3 %)ifc(Qjfc( m )fcfc(m)n 
(NiNiFNaNbNdNefVacVbeVdf Ei< m ,^ m ,i (M™ fc /s ) ~ 3 

(NaN b V ab ) 2 NiN d N jk N fn V d f EKm^m^i^^sl^Vlfc^tnt^mJ.V), + ) 



+ ^ (NiNtfNaNbNcdNefVacVbeVdf J2l<m,s^ m ,l(hk m kl s )- 3 

.; (NaN b V ab ) 3 NiN k N jm N hn Ei<m, S #m,;( fc i fc '" fc is)~ 3 fc(Qifc(;)jfc(m)fcfc(m)^fc(; S )mfc(i S )n 



{NiNipNaNbNniNefPacVbeVdf T,l< m ,s^ m A k lk m k ls 



+9< 



The coefficients t; can be read off directly from the expression in (1331) . Then, the 
t nl parameter can be written as: tnl = (1 + 6) t nl = 0- + + X2)t"nL' where 

is the part of tnl that depends only on the isotropic TS in (|2"4"|) and the isotropic 
PS in (|12p. The ratio again depends on the g^ parameter and on the configuration 
described by the ki momenta. 

5. A consistency relation for /nl and tnl 

Searching for relations between the NG parameters, we start exploring their 
isotropic pieces. Here we give an explicit proof that the isotropic pieces of the 



(33) 
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NG parameters, in presence of scalar and vector field perturbations, obey the SY 
consistency relational Using the expressions in (fl"2j) . (fT8|) . and (|24|) we find that: 



6 (to) _ N a N b N cd V ac V bd 
5 /nl (N a N b V ab f ' 



(too) _ N a N b N cd N ef V ac P be V df 



Now, we can write the scale invariant power spectra V ab as: 



V.,,, = D„uP ,<,, = ( J J ) P.,.. i3(>) 



where s = V+ /P54, is the vector to scalar PS ratio. Defining the vectors K a = D ab N b 

,-iso / / 6 f too ^ 
r NL/ \ 5 J NL J 



and T a = N ab K bl we evaluate the ratio A Nh = r^ s °/(f /^ s °) 2 : 



_ (N a K a )(D ab T a T b ) _ (D ab N a N b ){D ab T a T b ) 
NL (ifJ„) 2 (A^Ti,) 2 



(37) 



which implies: 



^>(^£) 2 . (38) 



The latter expression is a direct consequence of the Cauchy- Schwartz inequality 
for the vectors N a and T a with the positive definite metric D ab , and reflects the 
fact that the SY consistency relation, at least at tree level, applies for the isotropic 
pieces of the NG parameters even in the presence of vector fields. 

Now, we can try to establish a relation for the complete NG parameters /nl and 
tnl • The fact that the isotropic NG parameters obey the SY inequality implies that 
the total NG parameters obey a modified consistency relation. The modification to 
this relation will depend on several aspects, for instance, on the level of statistical 
anisotropy, the momenta configuration, etc. Employing the SY inequality for the 
isotropic parameters, we deduce that the total NG parameters obey the relation: 

(i + 6)(i + xi)V6. J 



TNL ^ 7T~ v-i , c n 2 7-/nl • (39) 

(l + X2)(l + £i) \5 / 

On the other hand, the relation obeyed by the NG parameters and r^ L is the 
following: 



(1 + 6) (G 

(1 + 6) 2 U - 



r NL > TT-7% ( f/nl ) • (40) 



All of these relations represent a constraint on the NG and statistical anisotropy 
parameters for models of inflation that involve vector fields. 
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6. Conclusions 

In this paper, we have studied under which conditions new kinds of the well known 
Sayama-Yamaguchi consistency relation are allowed when vector fields are included 
in the inflationary dynamics. In particular we have shown that the isotropic pieces 
of the non-gaussianity parameters obey the Suyama-Yamaguchi consistency relation 
t nl ^ (|/nl) 2 wnen its calculation is performed at tree level. 

We have also derived a modified consistency relation for the complete non- 
gaussianity parameters /nl and tnl and for a set of new parameters f^ L and r^j L 
that could be observationally relevant; these relations are given by the Eqs. (|39j) 
and PDI . Since the latter depend on the level of statistical anisotropy and the 
configuration of the wavevectors, the naive relations tnl (|/nl) 2 and r^ L 
(|/nl) 2 cou ld eventually be violated in some particular inflationary model that 
includes vector fields. On this basis, we expect that in a near future the derived 
relations could provide a criteria to discriminate among the different models for the 
generation of the primordial curvature perturbation when vector fields are involved. 
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